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Abstract 

We report on new results of the spin dependent structure functions g\ and 
h\ of the nucleon. An attempt is made to convert the moments, which is what 
one computes on the lattice, to quark distribution functions. 



1 Introduction 

In the past polarization data have often been the graveyard of fashionable models. 
Measurements of the polarized deep-inelastic structure functions of the nucleon over 
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the last decade have once again borne out this experience. While the naive quark 
parton model has been very successful in predicting the gross features of hadrons, it 
was a surprise to find that it fails to explain the spin properties of the nucleon. 

There is significant interest in an ab initio calculation of the nucleon structure 
functions now. The theoretical framework of such a calculation is the operator prod- 
uct expansion. While the Wilson coefficients can be computed perturbatively, the 
hard part of the calculation is the determination of the forward nucleon matrix el- 
ements of the operators. This is a non-perturbative problem, and the technique to 
solve it is lattice QCD. For recent work on the subject see A confrontation 
of the experimental data with the lattice predictions will be a crucial test of our 
understanding of the structure of the nucleon. 

At the twist two level the quark sector of the nucleon is completely specified by 
the spin averaged structure functions Fi(x, Q 2 ), F 2 (x, Q 2 ) and the polarized structure 
functions gi(x,Q 2 ), hi(x,Q 2 ). In this talk we will report on new results of the 
polarized structure functions. 

The talk is organized as follows. In sec. 2 we introduce the notation. In sec. 3 we 
give the results for the moments of g\ and h\. We furthermore give an extrapolation 
of the axial vector coupling of the nucleon, g^, to the continuum limit. The moments 
can be converted into real structure functions by an inverse Mellin transform. We 
will present first results on gi(x,Q 2 ) in sec. 4. Finally, in sec. 5 we conclude. 



2 Basics 

The structure functions g\ and h\ have simple parton model interpretations. The 
structure function g\ measures the quark helicity distribution in a longitudinally 
polarized nucleon, while h\ measures the probability of finding a quark in a spin 
eigenstate of the operator /175 in a transversely polarized nucleon. 

In the following we shall consider only non-singlet structure functions and distri- 
butions. Only these distribution functions are accessible in quenched lattice QCD. 

Let us first consider the structure function g%. If the sea is assumed to be flavor 
symmetric, the leading twist contribution can be written 

9i^Q 2 ) =4£ T ^ (*,%) Ag(/W), (1) 
z j Jx y \y H J 

where 

AqW (x, fi) = q\ f) (x, 12) - q{ f) (x, fi), (2) 

and qtf\ is the probability distribution of a quark of flavor / and spin parallel (anti- 
parallel) to the parent spin of the nucleon. The so-called splitting functions are 
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determined by the Wilson coefficients 

through an inverse Mellin transform. Similarly, the distribution functions can be 
derived from their moments, 

A (n) g (/) (/i) = f 1 dxx n Aq^(x,fj l ). (4) 
Jo 

We denote the lowest moment by 

Aq (f) = A(°y/). (5) 

According to the operator product expansion the moments are given by forward 
nucleon matrix elements of local operators. For A^ n 'q^' we have 

(£ ^ { *n}\Pi S> = ^AWg^K^ ■ • -P, n + ■■■- traces], (6) 

where 

/ i \ n 

lt-,n} = (2) ^Tfe D * ■ ■ ■ D ^ - traCeS - 

Here {• • •} means symmetrization. The lowest moment Aq^> measures the axial 
vector charge of the nucleon. 

Similar expressions can be derived for the structure function h\ 0. One simply 
has to replace A.qW(x,fj>) in eq. ([[]) by 

6qW(x, f i)=q ( P(x,ii)-q¥ ) (x, f Ji), (8) 

( f) 

where 9_j_(T) is the probability distribution of a quark of flavor / and spin ^±75 parallel 
(anti-parallel) to the spin of the nucleon, giving 

SWqV)( fi )= [ dxx n 8q U \x^). (9) 



We denote the lowest moment by 

5q U) = 5 (o) q (f). (10 ) 

The moments are given by the matrix elements 

2 

(P> s\°iZi-^}\P> $ = ^V^K^i - s^Pa)P2 ■ ■ -P^ + traces], (11) 

where 

(i \ n— ^ 
2) 9 (/) ^!75 ^- B^ n - traces. (12) 

The lattice operators, which are computed at the scale 1/a (a: lattice constant), 
must be renormalized and brought into the same scheme in which the Wilson coeffi- 
cients have been calculated. Generically we can write 

O i (ji) = Z ij O j (a), (13) 
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Figure 1: The quenched moments Au and Ad plotted as a function of a 2 . The lattice 
spacing is given in units of the string tension, K. The lattice data are denoted by •. 
the extrapolated values by O. The phenomeno logical values § are denoted by ^. 
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where the indices distinguish between the various lattice operators which can (and 
do) mix under renormalization. It is a task of its own to compute the renormaliza- 
tion constants accurately. For details of the calculation we refer the reader to the 
literature H and the talk of Schiller ||. 

The lowest moment 5q^' measures the tensor charge of the nucleon. In the non- 
relativistic quark model axial vector and tensor charges are equal, i.e. 5q^ = Aq^\ 
The structure function hi has opposite chiral properties to g\. It can be measured 
in polarized Drell-Yan processes, but not in deep-inelastic lepton-hadron scattering. 
For other ideas see ||. 

3 Moments 

The quantities to be computed on the lattice are the moments (0) and ([]). We have 
new results on Aq^ and gA, on A^q^ and on 8q^>. In the quenched approximation, 
which we are using, f = u,d, and we write and q^ = d. In the following we 

drop the argument /i from the quark distributions. 

Aq 

In lattice calculations one imagines the high frequency modes higher than the cut-off 
being integrated out. The logarithmically singular contributions of these modes are 
absorbed into the bare parameters, such as the coupling constant, while the power- 
behaved contributions are usually left unaccounted for. For Wilson fermions, which 
are widely used, the power corrections are of 0(a). By improving the action, these 
corrections can be reduced to 0(a 2 ) |J, thus giving results which are closer to the 
continuum limit a = 0. 

Removing 0(a) effects from the matrix elements requires improving the operators 
as well. For the axial vector current the improved operator, and its renormalization 
constant, have been derived recently J/J. We will make use of these results here. 

We have done calculations at two values of the coupling, (3 = 6.0 and 6.2, cor- 
responding to lattice spacings of a ~ 0.1 and 0.07 fm, respectively. The results for 
Au and Ad are plotted in Fig. [I]. The values given refer to the chiral limit, i.e. zero 
quark mass, which are obtained from the lattice data by a suitable extrapolation. 

As the remaining errors are of 0(a 2 ), we may fit the cut-off dependence by a 
formula of the form cq + C2fl 2 , and use this formula to extrapolate the result to the 
continuum (a = 0) limit. The outcome is shown by the solid lines. Clearly, we would 
have liked to have at least one more data point at another value of the coupling. 
We are working on that. We compare our results with the phenomeno logical valence 
quark distribution functions ||. The agreement is good, considering that the errors 
on the phenomeno logical values are of the order of 10%. 
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Figure 2: The axial vector coupling of the nucleon as a function of a 2 . The lattice 
spacing is given in units of the string tension, K. The lattice results are denoted by 
• . the extrapolated values by O. The experimental value is denoted by 
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9A 

A quantity, which is known very precisely experimentally, is the axial vector coupling 
of the nucleon, 

g A = Au - Ad = 1.26. (14) 

In Fig. ^] we show our results for g^, together with the extrapolation of the lattice 
data to the continuum limit. We find the continuum result to be in good agreement 
with the experimental value. 

This figure indicates quite clearly how important it is to correct for finite cut-off 
effects. Even after improving the action and the operators, cut-off corrections can 
be quite substantial still. We see that the lattice result increases by approximately 
20% going from our coarsest lattice at (3 = 6.0 to the continuum limit. 



Sq 

The moments Su, 5d have been computed using the improved action, but so far the 
operator has not been improved. We are currently working on this problem. For the 
renormalization constant we have taken the tadpole improved one-loop pertur- 



bative result JT0[. This means that the 0(a) corrections have not completely been 
removed in this opposed to the previous case. The result of the calculation is 

shown in Fig. [3]. This work adds another value of (3 to our previous calculation 



Before we discuss the result, let us find out what the tensor charge actually tells us 
about the structure of the nucleon. For a stationary nucleon the operator (12) differs 



from (|7|) by a factor of 70. In the non-relativistic limit fermions are in eigenstates of 
70, and so 5q and Aq are equal. By comparing 5q and Aq for a real nucleon, we can 
gain insight into how relativistic the constituents are. 

Comparing 5q and Aq in Figs. [I], H now, we see that they are equal within the 
error bars. This shows that a non-relativistic description of the spin structure of 
the nucleon is quite adequate. Does this mean that the quarks are in a relative s- 
wave and the missing spin is coming from the gluons and the sea quarks? Further 
investigations will have to show. 

Because the operator ([12]) with n — 1 is odd under charge conjugation, sea quarks 
do not contribute to 5q. This means that we might hope that the quenched calculation 
is giving an answer close to the true value. 



Another quantity which receives contributions from the valence quarks only is the 
second moment A^q. This moment is found from the n = 1 case of eqs. (j^j, (|7|). 
Again, the reason is that the operator is odd under charge conjugation. 
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Figure 3: The quenched moments 5u, 5d as a function of a 2 . The lattice spacing is 
given in units of the string tension, K. The lattice results are denoted by •. the 
extrapolated values by O. The numbers are renormalized at the scale /z 2 = 4 GeV 2 . 
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The special feature of this moment is that it is directly accessible experimentally, 
which allows us to test the valence quark distribution without further phenomeno- 
logical analysis. 



For the polarization asymmetry |L2| of n + minus ir inclusive cross sections one 
finds to lowest order in a, 



A7r+ _ 7r - = 4Au™\x) - Ad™\x) 
p 4u™ l (x) - d™ l (x) 



for a proton target, and 



= Au™\x) + Ad™\x) 

d M val (x) + d val (x) 

for a deuteron target. The fragmentation functions, as well as the sea quark con- 
tributions, drop out because of isospin invariance relating the various fragmentation 
functions with each other. 

The polarization asymmetries have been measured by the SMC-Collaboration. 
For the lowest moment of the valence quark distribution they found |H| Aw val = 
1.01 ± 0.19 ± 0.14 and AcP al = -0.57 ± 0.22 ± 0.11. The experimental errors are 
still a little too large to make a quantitative comparison with the lattice results. The 
SMC-Collaboration has recently extended their analysis to the second moment fL4 |. 
At fi 2 = 10 GeV 2 they obtain the result 

A (i) M vai = Q.169 ±0.018 ±0.012, (17) 
A (i) d vai = _ .055 ± 0.027 ±0.011. (18) 



A recent lattice calculation, renormalized at the same scale, gives (15 



A (1) m = 0.189 ±0.08, (19) 
A^d = -0.0455 ± 0.0032. (20) 

The lattice and experimental results agree within their respective errors. 



4 Distribution Functions 

The x-dependence of the structure functions carries valuable information about the 
dynamics of quarks and gluons which is not immediately available from the mo- 
ments. Furthermore, because of limited experimental data, moments are sometimes 
hard to compare with experiment. This is, in particular, the case for the higher 
moments. Theoretically, the structure functions can be obtained from the moments 
by an inverse Mellin transform. 

A first attempt of constructing nucleon structure functions from a few lower mo- 



ments was reported in [If:] for the unpolarized case. In this talk we shall consider a 
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Figure 4: The lattice results for the first three moments of Au(x) and Ad(x) at 
j3 = 6.0 and fi 2 = 4 GeV 2 . The curves are fits to the lattice data of the form (|2T|). 
The parameters of the fit are a = 0.04(6) and = 2.21(9). 
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Figure 5: The distribution x[Au(x) — Ad(x)], together with the phenomenological 
valence quark distribution. The solid line is the result of the fit, the dashed line is 
taken from ref. ||. 
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different method and apply it to g\. 

We restrict ourselves to the coarser lattice at (3 = 6.0. We have computed the 
three lowest moments. Calculation of a fourth moment is possible. The moments 
can be well described by the formula 

A (n) g = cr(/3)(n + l + a)- /3 . (21) 

This formula seems also to describe the moments of the unpolarized structure func- 
tions well. A fit is shown in Fig. |4j. An inverse Mellin transform of fl2l"D gives 

Aq{x) = cx a {-lnx) (3 - 1 . (22) 

In Fig. [| we compare the result for Au(x) — Ad(x) with the phenomenological distri- 
butions. The outcome is encouraging. We have hope that with one more moment and 
precise lattice data we are able to derive phenomenologically useful quark distribution 
functions. 

It must be said that one can only combine even and odd moments to make a 
single structure function if u and d sea quark contributions are assumed to be equal, 
as is commonly done M . However, this may not be a good approximation Jl7| . 



5 Conclusions 

Lattice calculations of nucleon structure functions have improved in many respects. 
The calculations are now done with improved actions and using improved operators, 
so as to reduce cut-off effects. Furthermore, the renormalization constants of the 
lattice operators, which are another source of errors, are gradually being computed 
non-perturbatively |7], 18]. On top of that, we have seen that it is important to do 



the calculation at several values of the coupling and do an extrapolation to a = 0. 
By adding one or two more data points, and with increased statistics, we will soon be 
ready to report reliable continuum results, at least in the quenched approximation. 

Where we can compare the lattice results with experiment or the phenomeno- 
logical analysis, we find good agreement. Our efforts over the last year have, in 
particular, paid off for g^, the axial vector coupling of the nucleon. Two years ago 
this quantity was considered a problem for quenched lattice QCD fl9fl . 
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